Existence of Orbital Order and its Fluctuation in Ba(Fei_2;Coa;)2As2 Revealed by 

X-ray Absorption Spectroscopy 



O 
(N 

o 

Q 
o 



o 
o 

Oh: 



o 



> 

(N 
(N 



X 



Y. K. Kini,i W. S. Jung,i G. R. Han/ K.-Y. Choi,^ K.-H. Kim,^ C.-C. Chen,3'4 t. P. Devereaux,^ A. Chainani,^ J. 

Miyawaki,5 Y. Takata,5'E| Y. Tanaka,^ M. Ouara,^ S. Shin,^'^ A. P. Singh/ J.-Y. Kim/ '[ll and C. Kim^'^ 

^Institute of Physics and Applied Physics, Yonsei University, Seoul 120-749, Korea 
^FPRD, Department of Physics and Astronomy, Seoul National University, Seoul 151-747, Korea 
^Stanford Institute for Materials and Energy Science, SLAC National Accelerator Laboratory, 
2575 Sand Hill Road, Menlo Park, California 94025, USA 
""'Department of Physics, Stanford University, Stanford, California 94305, USA 
^RIKEN SPring-8 Center, Sayo-cho, Hyogo, 679-5148, Japan 
^Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581, Japan 
^Pohang accelerator laboratory, Pohang University of Science and Technology, Pohang 790-784, Korea 

(Dated: December 13, 2011) 

We performed temperature dependent X-ray linear dichroism (XLD) experiments on an iron 
pnictide system, Ba(Fei_a;Co2,)2As2 with x—0.00 and 0.05, to experimentally verify existence of 
orbital ordering (00). We observed clear XLD in polarization dependent X-ray absorption spectra 
of Fe L edges. By exploiting the difference in the temperature dependent behaviors, we were able 
to separate GO and structure contributions to XLD. The observed 00 signal indicates different 
occupation numbers for dyz and dzx orbitals and supports existence of a ferro-OO. The results 
are also consistent with the theoretical prediction. Moreover, we find substantial OO signal above 
the transition temperature, which suggests that OO fluctuation exists well above the transition 
temperature. 

PACS numbers: 74.25. Jb,74.70.Xa,78.70.Dm 



Various experimental results have shown that there are 
anomalous in-plane anisotropic behaviors in iron pnic- 
tides well below the transition temperature [l|, A 
large anisotropy in spin flip energies Ja and Jb was 
obtained from inelastic neutron scattering experiments 
on CaFe2As2[l|. Quasi particle interference patterns 
in scanning tunneling microscopy result also suggest 
a large anisotropic behavior in electronic structure of 
CaFe2As2[l]. However, such large anisotropic proper- 
ties were not observed in other early experiments such as 
angle resolved photoemission spectroscopy (ARPES) 3- 
0], X-ray absorption spectroscopy (XAS)p| and optical 
measurements (IR) ff\ . 

This inconsistency comes from the existence of twin 
domains 0. However, it did not take long to devise a way 
to remove the twin domains Q and anisotropic behavior 
was indeed obtained in transport (il. [loj . ARPES [ill-[l3 
and opticalQ measurements. Observed anisotropy is 
not only quite large but also anomalous. Resistivity mea- 
surements reveal that the antifcrromagnetic (AFM) di- 
rection is more conductive than the ferromagnetic (FM) 
direction (lo|. which is opposite to the common expecta- 
tion that AF ordering suppresses the electrical conductiv- 
ity. This anomalous anisotropy cannot be explained by 
the structural anisotropy of the orthorhombic phase ei- 
ther. The difference in in-plane lattice parameters of the 
orthorhombic phase is less than 1 % 15| and such a small 
difference is not expected to explain the large anisotropy. 

An additional effect is needed to explain such 
anisotropy and it was subsequently proposed that orbital 
ordering (00) may play a key role 16l42l|. Observed 



anisotropy was between the AFM (often referred to as 
a— or X—) and FM (6— or y—) directions, which suggests 
that only dy^ and dzx orbitals are important among the 
five Fe 3d orbitals since other orbitals have R4 symme- 
try. In this respect, among the various possible orderings, 
ferro-OO which occurs through unequal occupation of dyz 
and dzx orbital states was propsed[l7. 3"21|. The pro- 



posed 00 is supposed to be responsible for the anisotropy 
in, for example, transport properties and also possibly 
drives the magnetic and structural transitionsjlq. It is 
therefore important to experimentally verify existence of 
such an 00. 

It was proposed that temperature dependent X-ray lin- 
ear dichroism (XLD) experiment may reveal the signa- 
ture of OOjlJl. The difference between XAS data taken 
with linear polarizations along the AFM and FM direc- 
tions (that is, XLD) is expected to have a specific spectral 
shape. To resolve the issue on the 00 in pnictides, we 
performed XLD experiments on detwinned mother (x—O) 
and underdoped (x=0.05) Ba(Fei_a;Coa;)2As2. Clear 
XLD signal was obtained from the results of both mother 
and underdoped compounds after the contribution from 
the structural anisotropy was removed. Both the 00 
and structural signals were found to exist well above the 
transition temperature, indicating strong fluctuation. 

Single crystals with two doping levels of x=0 and 0.05 
were synthesized by the self-flux method[24]. Both struc- 
ture (Ts) and magentic (Tat) transition temperatures of 
the mother compound were measured to be 135 K while 
underdoped sample showed Ts and Tn at 67 and 56 K, 
respectively (see the inset in Fig. 1(b)). Detwinning was 
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FIG. 1: (Color online) (a) Fe Lz,2 edge XAS spectra from 
mother compound BaFe2As2 with two different polarizations 
E II X (blue) and E || y (red) at 30K. The lower curve is the 
XLD spectrum multiplied by 10 (black). The inset shows the 
experimental geometry and axes information, (b) Fe L2 and 
L3 edge XAS and XLD spectra for an underdoped compound 
Ba(Fei_a;Coa;)2As2 taken at 40 K. The inset shows the tem- 
perature dependent resistivity. Resistivity was measured on 
a twinned crystal. 



achieved by applying mechanical stress on samples fll'|. 
XAS experiments were performed at beam line 2A of 
Pohang Light Source (PLS) and at 17U of SPring-8. 
Samples were cleaved in situ under a pressure better 
than 1.0 X 10~^ Torr and were transferred immediately 
into the measurement chamber at a pressure better than 
7.5 X 10"" Torr. Spectra were taken at various temper- 
atures. All the absorption spectra were recorded in the 
total electron yield mode and were normalized by the in- 
cident photon flux at a gold mesh. The direction of light 
polarization was controlled by rotating the sample. 

Figure 1 shows Fe L^^2 edge XAS spectra taken with 
polarizations along the FM (y— ) and AFM (x— ) direc- 
tions as well as XLD signal (difference) for (a) mother and 
(b) underdoped compounds. The experimental geometry 
as well as the definition of the axes are shown in the in- 
set of Fig. 1(a). Axes can be defined with respect to the 
direction of applied stress. To obtain the XLD signal, all 
XAS spectra were normalized using the following proce- 
dure. First, linear background was removed to make the 
tail part parallel to the energy axis. The spectra were 
then normalized by the area normalization method. Af- 
ter these steps, we could obtain the XLD spectrum by a 
simple subtraction of two normalized spectra. The result- 
ing spectra show small but clear XLD signal. With linear 
polarizations along the x— and y— axes, the only possi- 
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FIG. 2: (Color online) (a) Schematic of the predicted OO and 
structure XLD signal 21] and (b) their expected temperature 
dependence. Total XLD (solid black lines) and 00 contribu- 
tion (dashed red) at various temperatures for (c) mother and 
(d) underdoped compounds, respectively. Normalized XLD 
defined as {Iafm-Ifm)/Il3 at ojq as a function of tempera- 
ture for (e) mother and (f) underdoped compounds. Struc- 
ture (blue line) and OO (red line) contributions to XLD for 
(g) mother and (h) underdoped compounds. 



ble dichroism signal comes from the difference between 
dyz and dzx orbital occupations because others have the 
same parity with respect to both xz— and yz— planes. 
Therefore, the observed XLD implies that occupations of 
dyz and dzx orbital states are indeed different. 
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However, one has to be careful in interpreting the XLD 
signal in Fig. 1 as being due to 00 because it also con- 
tains contributions from the structural anisotropy. The 
system has an orthorhombic structure with the broken 
R4 symmetry below Tg, and thus we may expect ad- 
ditional XLD purely from symmetry reason without an 
00. Especially, the predicted structural contribution has 
a behavior opposite to that from 00 (see Fig. 2(a)), 
which may make the XLD signal small. 21| Therefore, 
we need to find a way to separate the two contributions 
in the XLD. We argue that the structure contribution has 
an abrupt change across the Tg and saturates right away 
unlike the 00 contribution which is aspected to slowly 
enhances (see Fig. 2(b))[23. The abrupt change for the 
structure contribution is justified from the fact that lat- 
tice constants also abruptly change at the T5 and do not 
show appreciable change well below the T5 even with an 
external stress [iJ, [2^. One can exploit such different 
temperature dependences of 00 and structure contribu- 
tions to separate out the 00 contribution as detailed 
below. 

We performed temperature dependent XLD experi- 
ment and the results are plotted in Figs. 2(c) and (d). 
The total XLD (black lines) has a clear temperature de- 
pendent behavior. Upon a closer inspection, the behav- 
ior is rather anomalous. Instead of monotonic decrease 
with temperature lowering, the XLD has a pronounced 
behavior near T5. In addition, the shape of XLD spec- 
trum changes dramatically around T5. Particularly, for 
the case of mother compound, the spectrum at 120 K 
not only shows more pronounced dichroism but also has 
opposite behavior near the edge compared to the 30 K 
spectrum. Such pronounced XLD signal around T5 is 
mostly from the structure contribution because 00 con- 
tribution has not yet set in. As the temperature is low- 
ered even further, 00 contribution finally comes in and 
somewhat cancels out the structure contribution, result- 
ing in a smaller overall XLD signal. 

However, actual structure contribution does not ap- 
pear abruptly at Ts as discussed above. To trace the 
detailed temperature dependence of structure contribu- 
tion, we plot the magnitude of XLD signal at ujq divided 
by L3 main peak intensity as a function of temperature 
in Figs. 2(e) and (f) for mother and underdoped com- 
pounds, respectively. In both plots, we see that XLD sig- 
nal increases within a certain temperature range around 
Ts (indicated by gray area) as the temperature is low- 
ered. This increase is due to appearance of the structure 
contribution. Note that the temperature at which the 
XLD signal starts to increase is higher than Ts- This 
implies that the structure transition also occurs from 
that temperature. The structure contribution finally sat- 
urates at a temperature below Ts where the XLD signal 
decrease again. Such broadening of the structure transi- 
tion temperature could be attributed to the effect from 
the applied external stress which was observed in recent 
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FIG. 3: (Color online) 00 contribution to the XLD signal for 
(a) mother compound and (b) underdoped compound. 
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We also noted differences in the temperature depen- 
dent XLD from the two compounds. In the underdoped 
case, XLD shows an down-turn behavior (between 120 
and 140 K) before it starts to increase while it is not 
seen in mother compound case. In addition, XLD in un- 
derdoped is significantly stronger than that in mother 
compound. These facts imply that 00 in underdoped 
compound starts to appear at a much higher tempera- 
ture than Ts and than in mother compound. XLD spec- 
trum of underdoped compound taken at 120 K also shows 
small dichroism signal that is consistent with the shape 
of 00 contribution. Detailed discussion on this doping 
dependence will be given later. 

With this difference in the temperature dependence, 
we successfully separated the 00 and structure contribu- 
tion from the total XLD. First, to extract structure con- 
tribution, we subtract the spectra taken above T5 (150 K 
for mother and 120 K for underdoped compounds) from 
that taken just below the transition temperature (120 K 
for mother and 80 K for underdoped compounds). We 
must mention that there is a finite change in the 00 
contribution between the two temperatures which must 
be considered in estimating the true structure contribu- 
tion. It is refiected in our data analysis but the details 
will not be discussed here. The resulting spectrum is the 
structure signal and is plotted in Figs. 2(g) and (h) (blue 
lines). Shapes for different compositions are similar and 
are consistent with the theoretical prediction [2lj. Then, 
the 00 contribution can be easily obtained by subtract- 
ing the structural contribution from the total XLD signal 
based on the decreasing behavior mentioned before. Re- 
sulting 00 contribution is also plotted in Figs. 2 (g) 
and (h) as red dashed lines. The behavior is again quite 
consistent with the predicted one schematically shown in 
Fig. 2(a) [111 . 00 contributions at other temperatures 
are overlaid in Figs. 2(c) and 2(d) (red dashed lines). 
Note that all the 00 spectra have a consistent shape 
and only the magnitude decreases as the temperature 
increases. Eventually, the 00 contribution disappears 
at around 170 K for mother and 140 K for underdoped 
compounds. However, 170 K spectrum from mother com- 
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FIG. 4: (Color online) Overlaid OO contributions at vari- 
ous temperatures for (a) mother and (b) underdoped com- 
pounds. Temperature dependence of the XLD magnitude for 
total (black), 00 (red) and structure (blue) for (c) mother 
and (d) underpdode compounds. 



pound shows a small dichroism signal but we attribute 
it to simal due to out-of-plane polarization as observed 
earlier [6j. 

We now analyze the 00 contribution in detail. We 
re-plot in Fig. 3 the data in Figs. 2(g) and 2(h). Both 
spectra show a dip feature at photon energy of 706.7 eV. 
On the other hand, a peak feature appears at different 
energies; a peak at 708 eV for the mother compound and 
a broad feature at 709.3 eV for the underdoped system. 
As we aligned the light polarization along the x— (AFM) 
and y— (FM) directions , spectrum taken with the x— 
(y— ) polarization detects dzx {dyz) orbital. Therefore, a 
dip in the data means dyz has higher density of states 
than dzx- This fact directly connects to unequal occu- 



pation numbers for dyz and dzx, i-e. dyz state is less 
occupied than dzx as theoretically predicted 17, 19l-21|. 
On the other hand, a peak feature on the higher energy 
side implies the opposite situation. However, we note 
that the area of the dip feature is larger than that of the 
peak feature for both of the compounds, which is again 
consistent with the theoretical prediction. Theoretically 
predicted XLD magnitude is around 3 % while our ex- 
perimental values are 4.25 % and 3.74 % for mother and 
underdoped compounds, respectively. 

Before we move onto discussing details of the temper- 
ature dependence, we wish to briefly touch upon possible 
contribution of magnetic origin to XLD. We believe mag- 
netic XLD is negligible considering the fact that there is 
an excellent agreement between theoretical predictions 
and experimental results of 00 signal. This view is also 
consistent with an earlier assertion that XLD signal from 
Ni^"*" is simply the property of wave functions of d-orbital 
states, which is drawn based on the comparison of exper- 
iments and atomic multiplet calculations (27|. 

For a detailed investigation of the temperature depen- 
dence, we plot 00 contribution again in Figs. 4(a) and 



4(b). A pronounced and monotonic temperature depen- 
dence is already seen. In Figs. 4(c) and 4(d), we plot the 
magnitudes of the structure and 00 contributions as a 
function of temperature along with the total XLD. Un- 
like the total XLD, partial contributions show monotonic 
behaviors for both compositions. A noticeable aspect of 
the data is that both structure and 00 contributions 
appear well above the transition temperature. For the 
structure contribution, its presence above the Tg could 
be attributed to an extrinsic effect due to the external 
stress for detwinning. 2^, 2^ 

Turning our attention to the 00 contribution, its ap- 
pearance above the transition temperatures indicates ex- 
istence of strong 00 fluctuation above the transition 
temperature. It is consistent with recent resistivity 
resul t llOl and anisotropic band shifts observed in ARPES 
data[l3|. In both results anisotropic behaviors appear 
above the transition temperatures around 150 K for un- 
derdoped samples. Such anisotropy above the transition 
temperatures was attributed to the effect of applied stress 
which makes spin fluctuation stabilized and pushes T/v to 
a higher temperature ^26']. Co-existence and a similar be- 
havior of 00 and spin fluctuations may suggest a close 
relation between the two. A possibility is that stabiliza- 
tion of 00 by the stress drives spin ordering through spin 
orbit coupling. 

On the other hand, the temperature range in which 
00 fluctuation exists increases with doping. It is about 
45 K for mother compound and 80 K for underdoped 
case. This indicates enhancement of the fluctuation ef- 
fect upon doping. The presence of 00 fluctuation as 
well as its enhancement upon doping could mean an im- 
portant role of 00 fluctuation for the superconductivity 
in iron pnictides. Indeed, it was recently proposed that 
00 fluctuation could act as a glue for the pairing in iron 
pnicitides. An important aspect of 00 fluctuation medi- 
ated superconductivity is that the superconducting gap 
has 8+4- symmetry 2^ 2^ . The gap symmetry, whether 

it is s+-|_ or s^ , is a crucial information in understanding 

the superconductivity. Our finding supports the discus- 
sion in a recent report 28[ and could explain why some 
systems seem to have s++ gap symmetry while others 
S-i , depending on which of 00 or spin fluctuation pre- 
vails. A comparative XLD studies of materials with s++ 
or S-|-_ gap symmetries could further resolve the issue. 
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